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Abstract 



The conjecture of a hidden Eio symmetry of M-theory is supported by 
the close connection between the dynamics of D — 11 supergravity near 
a spacelike singularity and a truncation of an one-dimensional cr-model 
with Eio symmetry where all representations beyond S 1 1/(10) level i = 3 
■ are omitted. If this conjecture is right, higher-level representations should 

' especially capture the dynamics of further M-theory degrees of freedom. 

Unfortunately, the level by level determination of Eio commutators which 
is necessary to extend the model to higher levels is both an involved 
and toilsome task that requires computer aid. In this work, some of the 
relevant problems are exposed and algorithmic methods are developed 
which simplify key steps in the determination of explicit Eio commutators 
(— j , at higher levels. As an application, we compute the commutator of the 

level-two six-form with itself. 



1 Introduction 

Since Witten's discovery |TZ] that the strong coupling limit of the five consis- 
tent ten dimensional superstring theories may be an eZewen-dimensional theory 
whose low energy regime is the unique known eleven-dimensional supergravity, 
there has been a lot of activity trying to get a better understanding of the 
nature of this strong coupling regime of string theory, as this is is a highly 
promising candidate for the complete unified description of the dynamics of all 
the fundamental forces. While a lot could be derived about this theory in the 
last decade, our picture is still quite incomplete, and this is mirrored by the 
preliminary incomplete name M-theory. 

As the study of symmetries - possibly hidden and showing up only in cer- 
tain limits of a model - has been extremely fruitful for the advancement of our 
understanding of fundamental physics, during the last century, often leading to 
a major breakthrough, any indication of a not yet well understood symmetry 



of M-theory should be considered as a potential important clue, and subjected 
to close investigation. While the existence of hidden exceptional symmetries of 
eleven-dimensional supergravity, especially the hidden global En-d symmetry 
that occurs in Kaluza-Klein reduction on a d < 8-dimensional torus, is well 
known (see e.g. [3] for an introduction), and speculations about a hidden global 
infinite- dimensional Kac-Moody symmetry also have been around for a long 
time ^21: the conjecture of a fundamental E\q symmetry of M-theory recently 
gained interest by the discovery that for the bosonic sector of eleven-dimensional 
supergravity in a setting close to a spacelike 'big bang' singularity (where spa- 
tially separated points become causally disconnected), the degrees of freedom 
of the metric show chaotic oscillations of BKL type [2] that can be described 
by relativistic billiards EI E] taking place in a region of the shape of the 
fundamental Weyl chamber of E\q. Intriguingly, a one-dimensional cr-model [7] 
constructed from a truncation of £10 to the first three levels of a SX(IO) decom- 
position JSj is able to capture and reproduce the dynamics of a correspondingly 
truncated expansion in terms of gradients of D = 11 supergravity in the pre- 
viously described limit. While this already is a highly nontrivial result, one 
may still be sceptical whether this correspondence may be extended to higher 
levels. At least, E w seems to contain all the fields necessary to represent higher 
gradients of the 3-form and 6-form potential as well as the 'dual' graviton. Be- 
sides those, there are many other representations that are conjectured to be 
linked to other dynamical (potentially even yet unknown) degrees of freedom of 
M-theory. First promising indications that such a description of M-theory via 
an E10 er-model may indeed work are given by the highly nontrivial observation 
that it may be reconciled (in a truncation to low levels) with both (massive) 
II A and IIB supergravity p3 > as well as the possible correspondence of 
S'L(IO) singlets at quite high levels to R 4 , R 7 ,. . . corrections of M-theory 

Clearly, it is essential to perform an extension of the construction beyond 
level three. As we will see, this is bound to require some quite involved algo- 
rithms, and hence this line of research may well benefit from scientific exchange 
with computer scientists specialized in symbolic algebra. Therefore, the present 
article is somewhat more verbose than one targeted exclusively at a string theory 
audience. While it is not expected of readers with a specialization in symbolic 
algebra to comprehend all of the physics upon first reading, at least the relation 
between the relevant concepts should become clear. While we will be concerned 
exclusively with the determination of commutators for the Borel subalgebra of 
E10 in this work, details on the construction of the model up to level three can 
be found in [Sj. 

2 Conventions 

In the following, we will make extensive use of irreducible representations of the 
A n = SL(n + 1) Lie algebra that are obtained by applying mixed-symmetry 
Young projectors on tensor powers of the fundamental vector representation. 
We use the usual convention of first symmetrizing over rows, then symmetrizing 



2 



over columns of a Young tableau, i.e. 
P 



B\A abc = const. • ({A abc + A bac ) - (a <-> c)) 



(1) 



It is furthermore convenient to lay out all the indices of a tensor carrying 
a single irreducible representation in a form that directly mirrors the shape of 
the corresponding symmetrizer's tableau. However, as we will mostly have to 
deal with short symmetrizations and long anti-symmetrizations in the analysis 
of £10, we will rotate our tableaux by a counterclockwise quarter-turn. Thus, 
from now on, we write e.g. 



c\d\e\f\g\h\i\ 



instead of ES 



b 


a \ 


c 




d_ 

e 




/ 




<L 
h_ 




j_ 





(2) 



for the 'dual graviton' generator at level three in the Ag decomposition of £10 • 
Despite this notation, which is rotated solely for obvious typographic reasons, 
we will still speak of the symmetrization of rows and anti-symmetrization of 
columns/ 

We will also extensively apply symmetrization and anti-symmetrization op- 
erations to groups of indices. While this is usually denoted with brackets or 
parentheses around groups of indices, such as in 1 



b' c' d' 



-^-a[bcd]e ■— ^ab'c'd'eOf, c d 

this ad-hoc notation has two important drawbacks: 



(3) 



• It is highly unpractical for more complicated situations with interlocking 
symmetrizations. 

• Positional information is used to convey information about symmetriza- 
tions. 

In particular, there is no reasonable way to translate an expression like 



y ra'd'g' cb' f' ce'e' 

X a 'b'c'd'e'f'g'd adg f> bS <> ce 



(4) 



to full bracketed-indices symmetrization notation. Such quantities arc, how- 
ever, bound to appear in the analysis of higher levels of hyperbolic Kac-Moody 
symmetries. 

1 As usual, we define anti-symmetrizers as projectors, e.g. 
8f e c f = i (<55<5e<5/ ± permutations) . 
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The second problem mentioned has an important conceptual aspect: there 
are two different ways to apply an (anti-)symmetrization to a given mixed- 
symmetry tensor. One has to discern between using the antisymmetry in indices 
of the same column to exchange boxes, i.e. 





* 




* 




X 


* 


*1*2 * jr 


* 


*2l*ll * 



(5) 



and applying a symmetrization to index names, as in: 



X 



a 




b 


c \d e 



V 

b 



,b' -^a.b 



a 




b 


c \d e 



-X 



[6] 




a 


c \d e 



(6) 



which would correspond to 'placing brackets around [ab]\ 

One can make the idea rigorous that the first operation acts on the right, 
while the second one acts on the left of the operation of naming tensor indices, 
hence we will denote the first kind of symmetrization as a box shifting sym- 
metrization, and the second one as an index relabeling symmetrization. Note 
that box shifting symmctrizations are conceptually tied to individual columns of 
a single (mixed-symmetry) tensor, while index relabeling symmetries are more 
naturally associated to entire sums, not individual summands, as they amount 
to introducing an extra S y factor. A large part of the algebraic complexity 
of the higher levels is related to the fundamental difference between these two 
concepts. 

We clearly need a new notation for index relabeling symmctrizations that 
overcomes the forementioned problems. We will use the convention to denote 
symmetrization in indices i,j,k,...,n by appending a formal (pseudo-)factor 
a(ijk . . . n) to the right of the term, and likewise anti-symmctrization in these 
indices by a formal factor a(ijk . . .n). This definition implicitly contains the 
rule to apply such factors one by one from the left to the right to a term, re- 
placing a(ijk . . . n) by , where i',j', k', . . . ,n' are fresh indices, and 
substituting free indices i i' : j \— > j', etc. in all factors to the left of it. 
If it helps to make a point, we will occasionally include trivial one-index sym- 
metrization factors like a(g) in our formulae. 

To give an example, applying a specific Young projector may now be written 

as: 

a 



d 



(term) 



- (term) 
o 



cr(ab) ■ a(bcd). 



(7) 



The numerical factor cp, 8/5 in our example, that makes a sequence of a tableau 
symmetrization followed by a tableau anti-symmetrization an idempotcnt pro- 
jector, is given by 



c P = 



_ (Ilrerows^ 7 *)-) ' (Ilcecolumns^ )-) 



n 



h£ hooks 



(8) 
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(For every box in the tableau, there is exactly one hook, which consists of this 
box and all that lie either in the same row and to the right of it, or in the same 
column and further down.) 

Note that with this notation, it is very easy to make proper use of both box 
shifting and index relabeling symmetries. For example, we have: 



X 

= -X 

= +x 
= +x 



a 




b 


c\d\e 


a 




b 


d\c\e 


c 




b 


d\a\e 


c 




a 


b | d\ e 



a(a, b, c) 
a(a, 6, c) 
a(a, 6, c) 
a(a, b, c) 



(box shifting) 
(index relabeling) 
(2x box shifting) 



(9) 



This raises the question how to define normal forms for such expressions. 
This is an involved issue that has many technical and algorithmic aspects and 
hence is discussed in detail in the appendix in section |A. 1.21 In brief, the idea 
behind the scheme used in this work is to treat all extra (non-overlapping) index 
symmetrizations on a tensor as being ordered lexicographically, then using them 
one after another to put as many lexically small tensor indices into columns to 
the far right as possible. This means in particular (as is shown with an example 
in the appendix) that submitting a mixed-symmetry tensor in normal form to 
its associated Young symmetrization will in general produce a sum of more 
than one tensor in normal form. This is mostly due to our ignorance of more 
complicated relations between mixed symmetry representations with permuted 
indices, which are most easily expressed in terms of Garnir symmetries |11| - 
these will be explained in some more detail at the end of section EP1 

Defining a normal form that takes into account all three types of column 
antisymmetrization, index relabeling, and Garnir symmetries, and can be com- 
puted with sufficient efficiency for an application to E\q is not as inextricable 
as it may seem at first, but will not be discussed here, even if using such a 
normal form that differs from the one used here will eventually turn out to be 
inevitable when going to high levels. Note that if one would not have to deal 
with index relabeling symmetries, it would be possible to define normal forms 
of filled tableaux just such that in a tensor in normal form, all indices appear 
in lexically ascending order in every column as well as in every row. Any tensor 
can be transformed into a linear combination of tensors of that form by repeated 
application of Garnir as well as column antisymmetries. It is easy to see that 
this prescription has to be refined further to get unique normal forms in the 
presence of index relabeling symmetries, as is demonstrated by the following 
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example: 



X 



[61 




a 


c | 



cr(6c) 



c 




a 





a(bc) 



(10) 



One can imagine that the situation will get far more involved in situations 
where e.g. index relabeling symmetries touch indices that are not located at the 
boundary of the tableau. 

A further issue concerns the use of epsilon tensors: in Young tableau index 
notation, an element of the Eiq algebra which appears in the n-fold iterated 
commutator [[. . . [Ef c , E( e / ], E{ h % E{ kl ] . . .] will carry 3n indices. Typically, 

V * ' 

n 

some of those will form columns of maximal length (here 10) that could be 
extracted as extra factors ( ahc - -i . For the sake of a homogeneous presentation, 
we will not do so, and instead keep all columns of maximal length explicit. 

3 The £"io Borel subalgebra from the A§ perspec- 
tive 

We start this section by reviewing a few ideas and facts from Cartan-Dynkin 
theory and its application to E\q. The roots of the E\q algebra lie on the 
unique even unimodular Lorentzian lattice 7/9,1. Taking the ten-dimensional 
scalar product to be 



(x,x) 



£4 

this lattice is spawned by these ten vectors of squared length two (cf. p^) 



(11) 



rg 



1 
2 




-1 

1 













1 1 

2 2 















-1 





(12) 



The Dynkin diagram encoding the geometrical structure of this lattice is 
given in figure ^ 



3.1 Kac-Moody Algebras from Cartan matrices 

One key property of Lie algebras is the geometry of their root lattice. Choosing a 
hyperplane through the origin that does not contain any further lattice point in 
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123456789 

Figure 1: The Dynkin diagram of Ag extended to E w 



order to separate lattice vectors into a 'positive' and 'negative' half and denoting 
positive lattice vectors that cannot be written as the sum of other positive lattice 
vectors as simple roots, this information can be concisely encoded in the Cartan 
matrix Aij , which contains the scaled (not necessarily euclidean) scalar products 
of the simple roots: 

A« = 2^ (13) 

The off-diagonal entries of a Cartan matrix are nonpositive integers, due to 
the indecomposability of the simple roots. Vice versa, it is possible to construct 
the Lie algebra qa from a generalized Cartan matrix: one first defines the 
Chevalley generators ej,/i,/ij, with commutators satisfying 

[hi,ej] = AijCj [hi, fj] = —Aijfj . . 

[e»,/j] = Sijhi [hi,hj] = 0. 

The algebra qa is then obtained first forming the free Lie algebra, which 
is spawned by all the commutators that can be built recursively from these 
generators, modulo the Jacobi identity 

[[X, Y],Z} + [[Y, Z],X} + \[Z, X],Y} = 0. (15) 

This free Lie algebra then has to be divided by the Serre relations: 

■ ■ ■ , [Si, Bj] 

S v ' 

[fi i \fii •■•;[/»> fj] 
v v ' 

For finite-dimensional Lie algebras, all roots have squared length A 2 = 2. In 
the hyperbolic case - where we will take E w as our primary example - the metric 
on the root space is indefinite, and while there is just one generator associated 
to every root with A 2 = 2 (as in the finite-dimensional case), the multiplicity 
of roots inside the 'light cone', with A 2 < 0, increases quite dramatically with 
increasing — A 2 . While the Jacobi and Serre relations determine the structure of 
the Lie algebra, it is a difficult and basically still unsolved problem to effectively 
identify those multiple commutators of Chevalley generators Ei that vanish due 
to the Serre relations, which may be revealed only after multiple applications of 
the Jacobi identity. 



•]] = 
.]] = 0. 



(16) 
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While one may construct a Lie algebra basis as well as determine the struc- 
ture constants directly from the simple roots, it is frequently more convenient 
and more useful to identify a simple subalgebra that has a simple geometrical 
description by omitting one (or more than one) simple root, and slice the alge- 
bra into irreducible representations of this sub-algebra, graded by the number 
of occurrences of the operator Ej or Fj corresponding to the omitted root rj . 
This leads to the presentation of commutators in the language of tensors of the 
sub-algebra. An example of this technique, the decomposition of Eg in terms of 
its SL(8) = A 7 sub-algebra, can be found in appendix B of 15 . 

In this work, we limit our analysis to E\q Borel subalgebra commutation 
relations that arc schematically of the form 

Ee =p+q = y^jEj-p, Ej— q ]. (17) 

Most of the technical complication comes from the demand to work not with 
explicit roots, but abstract tensor equations, and the combinatorical infeasibility 
to use naive (anti-)symmetrization by generating all permutations if more than 
eight indices have to be (anti)-symmetrized. Hence, the question arises how to 
organize calculations in such a way that maximal use can be made of overlaps 
with pre-existing symmetries. 

3.2 The A 9 decomposition of E 10 at level £ < 3 

The determination of commutators for the E w algebra not only gets more and 
more voluminous with increasing Ag level, one also has to use an increasing 
number of different algebraic tricks with increasing level. Also, some steps per- 
formed for the first few levels of the construction can be interpreted in multiple - 
sometimes unexpected - ways, and often, only one interpretation can be suitably 
generalized to be of use for higher levels as well. 

Hence, it makes sense to start the discussion of higher levels with a brief 
review of the derivation of commutators up to level £ — 3, where we mark those 
steps in the construction that require closer inspection, resp. modification, when 
going to higher levels with (*„). A more detailed discussion of the nature of 
these steps will then be given in section I5~4l 

One important tool in the derivation of commutation relations is the repre- 
sentation table in |16| . whose first few entries we quote for convenience: 



Table 1: Ag representations in E w up to level I — 6 



t 


V 


771 


A^ 


dim K(A) 


mult (A) 


/J 


1 


(001000000) 





























2 


120 


1 


1 


2 


(000001000) 


1 


2 


3 


2 


1 














2 


210 


1 


1 


3 


(100000010) 


1 


3 


5 


4 


3 


2 


1 








2 


440 


1 


1 




(000000001) 


2 


4 


6 


5 


4 


3 


2 


1 








10 


8 





4 


(001000001) 


2 


4 


6 


5 


4 


3 


2 


1 





2 


1155 


1 


1 




(200000000) 


1 


4 


7 


G 


■5 


4 


3 


2 


1 


2 


55 


1 


1 



8 





(010000000) 


2 


4 


7 


6 


5 


4 


3 


2 


1 





45 


8 





5 


(000001001) 


3 


6 


9 


7 


5 


3 


2 


1 





2 


1848 


1 


1 




(100100000) 


2 


5 


8 


6 


5 


4 


3 


2 


1 


2 


1848 


1 


1 




(000010000) 


3 


6 


9 


7 


5 


4 


3 


2 


1 





252 


8 





6 


(100000011) 


3 


7 


11 


9 


7 


5 


3 


1 





2 


3200 


1 


1 




(000000002) 


4 


8 


12 


10 


8 


6 


4 


2 








55 


8 







(010001000) 


3 


6 


10 


8 


6 


4 


3 


2 


1 


2 


8250 


1 


1 




(100000100) 


3 


7 


11 


9 


7 


5 


3 


2 


1 





1155 


8 


1 




(000000010) 


4 


8 


12 


10 


8 


6 


4 


2 


1 


-2 


45 


44 


1 



At level one, there is the single three-form corresponding to the root tq. As 
there is only a single six-form at level two, the only sensible level two commutator 
that can be written down (*i) is 



[Ef c ,Ef ef ] = aEf cdef . 



(18) 



As the normalization of E% can be chosen arbitrarily, we set the coefficient 
a=l. 

The level three, we get a single tensor of mixed symmetry. Again, the most 
general ansatz for the level three commutator is 



[Ef cdef ,Ef hi ] 



c I d | e | f | g | h \ 



b|c|d|e|/|/i|z| 



■a(abcdef)-a(ghi) 



(19) 

as there are only two different ways how to distribute indices from two columns 
of lengths (6,3) to two columns of lengths (8,1). If we make use of the identity 



P ( \a\b\c\d\e\f\g\h\i\ ) e 
a 



c\d\e\f\g\h\i\ 







6E 



c I d | e | f | g | h \ 



' + 33 



b\c\d\e\f\h\i\ 



■ a(abcdef) ■ a(ghi) = 
(20) 



and the freedom in the normalization of £3, this can be reduced to 



abedef j^Qhi 



b\c\d\e\jj\h\i\ 



a(abcdef) ■ a(ghi). 



(21) 



In order to obtain an expression for E3 without the extra index relabeling 
symmetrizations in terms of commutators [E\ , E 2 ] , we first apply an index rela- 



beling symmetrization (* 2 ) enforcing the column anti-symmetry of 



b\c\d\e\f\h\z\ 
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to the entire expression: 



[E2 bcde f , Ef hi ] ■ a(g) ■ a(abcdefhi) 



b\c\d\e\f\h\i\ 



2B 



\d\e\f\g\h\ 



a{g) ■ a(abcdefhi). 



(22) 

From this point on, we will only work with terms with the proper column 
symmetry a(g) ■ a(abcdefhi). 

Besides the desired term, we get an extra contribution where the tableau 
structure does not match the index symmetrizations. This can be eliminated 
by making use of the Jacobi identity (#3) in the form 



[[Ef c , £f f ],El hi ] + {[E™ 1 ,E B x m ],E1 bc ] + [[£?'", Ef c ], E™ 1 ] = 0. (23) 

which can be reduced to an equation between £3 tensors by use of the defini- 
tions (|21I32|) 



^def rpghii ™6c 



ghi j-iabci jpdef-i 



\d\e\f\g\h\ 



- B 



b\c\e\f\g\h\i \ 



b\c\d\e\f\h\i\ 



(24) 



■a(abc) ■ a(def) ■ a(ghi) = 0, 
which then has to be brought to a(g) ■ a(abcdefhi) symmetry: 



8B 



c\d\e\f\g\h\i\ 



b\c\d\e\f\h\i\ 



a{g) ■ a(abcdefhi) = (25) 



V 



This can then be used to reduce l)22|) to 



E abcdef^ E ghi^ . ^ . a ( abcde f M j = 1£ 



b|c|d|e|/|/i|z| 



(26) 



which is of the desired form. 



3.3 E w at A, level 



When going to level £ — 4, a few new features are encountered. First and fore- 
most, the commutator [E 3 ,Ei] will now contain two different representations. 
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The most general ansatz which corresponds to (I19f) for level I — 3 now is 



c|rf|e|,f|fl|/*Rl gjkl-i 





j 


k 


I 




aiE 


a 


b 


c 


d\ 


e\f\g\h 




4 














a 


3 


k 








a 3 E 


b 


c 


d 


e 


f\a\h\i 


















a 


b 


J 








a 5 E 


c 


d 


e 


f 


g\h\i\k 


















a 


b 


c 








ot-jE 


d 


e 


f 


f) 


h\i\j\k 


1 1 


4 

J_ 




PiE 


a 


b 


c 


d 


e\.f\g\h 


i\l\ 


4 

b_ 

j 






a 


c 


d 


e 


f\g\h\i 


fc I 



a 2 E l 



a^E l 



b 


J 


fc 




- a 


c 


d 


e\f\g\h\i\l 










b 


c 


3 




r- a 


d 


e 


f\g\h\i\k\l 



b 


c 


d 




a 


e 


f 


g\h\i\j\k\l 



(27) 



c|d|e|/|ff|/i|i|fc|Z| 



d I e | / | g | h \ i \ j \ k \l\ 



d | e | / 1 ff | h | i | j | fc | I | 



■ a(a) ■ a(bcdefghil) ■ a(jkl) 



J 



taking into account all possible ways to re-distribute indices from the three 
columns (a;bcdefghi;jkl) to two columns of length (3,9), resp. three columns 
of length (1, 1, 10). Clearly, we have (3i = 0, due to the anti-symmetry in jk. 
This approach, however, does not yet fully appreciate the mixed symmetry of 



the tensors on the right hand side. Thinking of E/^~ 
linear combination of basis tensors that diagonalize the regular part of the cor- 
responding Young projector, we can use the fact that sequential application of 
any mixed-symmetry projection (i.e. any arbitrary distribution of the indices 
to symmetrize and anti-symmetrize over) corresponding to a tensor of this par- 
ticular form to any such tensor, potentially with a different index distribution, 
followed by Young projections to the symmetries of the tensors on the left hand 
side, yields either zero or (up to an integer multiple that can be absorbed in the 
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normalization) just one specific combination: 



54i 





a 




a 


.7 


k 




J.7IAI/L 


b 


c|d|e|/|<?|ft|iL 


b 


c 


d 


e\f\q\h\i\l 



<E 



a 


b 


c 




d 


e 


f 


a\h\i\j\k\l\ 



81 j \i\k\l\Ab\c\d\e\f\g\h\i\ 



a 


.7 


k 




a 


.7 


k 




b 


c 


d 


e\f\g\h\i\l\„ 


b 


c 


d 





a 


b 


c 




a 


b 


.7 




d 


e 


f 


a\h\i\j\k\l\ , _ 


c 


d 


e 


/MfcUlfcU 



o 


.7 


k 




6 


c 


d 


e\. 




i\i\l\ 



b 


c 


d 




& 


c 


,7 




a 


e 


f 


fl|/lUI.7|fcU| 


fl 


d 


e 


/Ifll/iUlfcU 



=: Zi 



Likewise for the symmetric rank-2 contribution: 



b 


,7 






a 


c 


d 


el.fl.Ql/iUUI 



(28) 



_ 108i JIEJ]p[I 

si jum 



■Irflel/lffl/iUI 



c|d|e|/| g |fe|i|A:r 



d\e\f\ g \h\i\j\kW 



c\d\e\f\g\h\i\ 



c\d\e\f\ g \h\i\k\ 



= -m 



d\e\f\q\h\i\j\kTT 



c|d|e|,f|q|/i|TTfcH1 
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Hence, if we start with 



(29) 



(30) 



we may in principle be able to derive the form of the [E 2 ,E 2 ] commutator 
by making use of the Jacobi identity in the schematic form [[E 2 , ^i*^]) E^] + 
[[E[ i} , E[ U) ], E 2 ] + [[e[ U \ E 2 ], E[ i} ] = 0. Before we discuss this somewhat subtle 
issue further, we give the result. The only relevant nonvanishing commutator 
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is: 
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■a(abcdefghi) ■ a(jkl) 



3.4 On the role of the Jacobi identity 

As the structure of the Lie algebra is determined fully by the Serre relations 
and the Jacobi identity, and both enter (indirectly) into the determination of 
root multiplicities via the Peterson formula, one may ask the question whether 
the problem that the Serre relations and Jacobi identity are interlocked in an 
unwieldy way may be circumnavigated by trying to trade the Serre relations for 
the known root multiplicities in the level by level determination of the structure 
constants. In other words, taking a commutator definition of the form of (|30|) . 
what information can be extracted from the Jacobi identity at the given level? 
In particular, is it possible to derive further relations of the form l|25[) that can 
be used to simplify the commutator? 

A head-on approach to this problem at level £ is to systematically generate all 
possible commutator structures of i tensors E\ and all linear relations between 
three of them determined by the Jacobi identity. Taking all possible index dis- 
tributions into account, and re-symmetrizing to a given desired Young tableau 
column antisymmetry, this is bound to produce all the relevant relations that 
can be obtained in such a way. Clearly, this procedure soon becomes infeasible 
even on a powerful computer due to combinatorical explosion, but it is neverthe- 
less important to try this for the following reasons: First of all, there are many 
opportunities where even a small mistake in the implementation of an extensive 
algebraic algorithm can lead to wrong results that are hard to detect with the 
naked eye, do not occur for many simple test cases, but lead to inconsistencies in 
large calculations. An extensive analysis like the one described is one of the few 
but crucially important 'smoke tests' available that help to discover program 
bugs. Unfortunately, it also can only give binary information: either simplifi- 
cation of a considerable number of linear equations (typically > 200 equations 
containing in total 1000 summands with twelve-index tensors at level four) gives 
nonsensical relations like a vanishing commutator that is known to be nonzero, 
or a set of equations that seem to make sense. In the former case, the so far only 
way to find a bug seems to be to print out a detailed trace of all term transfor- 
mations in the program and check this manually - which easily becomes highly 
frustrating sysiphus labour in the long run. The second reason is that this is the 
conceptually simplest approach to the systematic determination of commutators 
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of the form [E mi Eg^ m ] for m > 1. While one might initially try to start with 
a definition for [E\,Ei-i] and use the Jacobi identity to derive the commuta- 
tor [E 2 ,Ee^2] by splitting £^_i to a term of the schematic form 
and then continue with this procedure to successively determine [E3, E1-3], etc., 
this approach is complicated considerably by the fact that a direct application 
requires using a definition of the schematic form Ei_ m = ^2[Ei, Ei- m -x], while 
one starts with an equation of the form [Ex, £!f_ m _i] = ^ E^ m , with tensors 
belonging to a variety of irreducible representations on the right hand side, that 
first has to be solved. Developing the relevant algorithmic techniques to effi- 
ciently automatize this is simplified greatly by knowing the right answer from 
the start from a conceptually simpler and hence more robust approach that may 
well be not efficient enough to be used for high levels. 

Efficient algorithms to generate a complete set of level I double-commutator 
identities from the Jacobi identity are described in appendices IA.2IA.3| These 
furthermore have to be instantiated with all possible distributions of the indices 
abedef ghi ... to the fundamental E\ tensors, modulo permutation of groups 
of three indices, e.g. (abc) (def), and then re-symmetrized in all possible 
different ways to the desired symmetry. Special care has to be taken with 
substitutions of commutator definitions: it is easy to overlook that an equation 
of the form [E t , E m ] ■ a{. . .) ■ a{. . .) — ^(tensors) • a{. . .) ■ a(. . .) may only be 
substituted into a term containing [E e , E m ] ■ a{. . .) ■ a{. . .) after suitable index 
re-labeling if the distribution of indices in the (anti-)symmetrizers a(. . .) relative 
to the distribution of those indices to the columns of Ee,E m can be made to 
agree. Efficiently finding suitable index relabelings so that a substitution can 
be applied is a further nontrivial task. 

Looking at the step from level £ = 3 to £ = 4, the major points that are not 
yet visible at level t = 3 and require generalization for an algorithmic approach 
are the following: 

(*i): For a simplified uniform systematic algorithmic treatment, all commu- 
tator definitions have to be given in a form such that both left hand side and 
right hand side make the column anti-symmetries of the tensors in the com- 
mutator explicit. Hence, even the [E\,Ei] commutator should be given in the 
form 

[Ef c , Ef ef ] ■ a{abc) ■ a(def) = aEf cdef ■ a(abc) ■ a(def), (32) 

whose anti-symmetrization redundancy might look slightly nonsensical, but ac- 
tually helps greatly to simplify many algorithms. 

(* 2 ): One should use a sequence of Young tableau projections instead of just 
the corresponding column anti-symmetrizations to reduce commutator defini- 
tions. 

(*3): At least at level t = 4, trying to use the Jacobi identity to derive further 
relations between 'wrongly symmetrized' tensors of mixed symmetry seems to 
be a red herring. The Jacobi identity does not produce any new relations of 
the strict form ^(tensors) = 0. Even at level £ = 3, it is perhaps more helpful 
to think of the corresponding identity as a consequence of a Garnir symmetry 
and not the Jacobi identity: for every pair of horizontally adjacent boxes in a 
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tableau, there is a linear relation between all tableaus with fillings that can be 
obtained by permuting the indices in those two boxes, plus all below the left 
one, plus all above the right one, in all possible ways. The coefficients in that 
linear relation are ±1, and determined by the constraint that the relative sign 
between terms that differ by an exchange of adjacent boxes is —1. Thus, for the 



two horizontally adjacent boxes in the tensor E 
relation: 



c I d | e | / 1 g | h \ z 
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-E 



\d\ 



a\d\ e\ f \ g\h\ 



-E 



\f\a\h\ 



we get the 



(33) 



= 



Anti-symmetrizing with a(abcdefhi) then gives l|25|) . Hence, if the notion 
of a tensor normal form already respected Garnir symmetries, the reduction 
from l|22l) to (|26[) would be automatic. 



3.5 Beyond level four 

There are a few obvious issues that have to be kept in mind when going from 
level four to level five. First of all, one no longer has one single commutator of 
the form [E4,Ei], but two different level four representations whose commuta- 
tors with Ei have to be taken. Furthermore, one no longer can just normalize 
-E5 tensors such that they appear with coefficient 1 in such a defining equation. 
The reason is that the 100100000 representation at I = 5 may a priori appear in 
the commutator of both £ = 4 representations with E\ . One might still argue 
that it will presumably at least appear in the commutator of the 200000000 rep- 
resentation with Ei, as this commutator cannot produce the only other £ = 5 
representation. If it were zero, the generators belonging to 200000000 would 
have to commute not only with E\, but (via the Jacobi identity) also with 
E2, and (inductively) with all generators at higher level. Hence, the coefficient 
could be set to one in this particular commutator. Such reasoning cannot be 
carried any further, however, as the 010001000 representation at level £ = 6 is 
in the tensor product of both £ = 5 representations with the level-one repre- 
sentation, which both carry further tensors. At level £ — 7, we encounter outer 
multiplicities /1 > 1 for the first time. 



4 Conclusion 

Lacking a powerful theory of the structure of hyperbolic Kac-Moody algebras, 
their possible relevance in gravitational physics - especially the conjectured role 
of E w for M-theory - probably justifies going to great lengths to obtain detailed 
information about even only the first few elements of a finite truncation. This 
is virtually impossible to achieve to satisfactory depth via manual calculation, 
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both due to the sheer number of steps and due to the inherent ineptness of 
humans to execute them in such large numbers without making mistakes. Hence, 
this task is bound to require computer aid. Unfortunately, the algorithmic 
side of this problem turns out as being of unusually high complexity, hence 
despite considerable technical effort, only very little could be achieved in this 
work in terms of new data about E 10 structure constants. The most important 
obstacle to further progress seems to be the unavailability of powerful tests that 
simplify the discovery and analysis of flaws and errors in algorithms and their 
implementations: while it is essential to try out various approaches to individual 
sub-problems (and, as it turns out, frequently change them when they have 
to be embedded into a larger context), all that can be done to validate an 
algorithm (not to speak of verification) is to use it to automatically generate all 
conceivable relations from a certain not too small set, semi- automatically check 
for inconsistencies, and - should they arise - eventually redo dozens of lengthy 
calculations by hand that involve tensors with more than ten indices, where 
one has to especially take great care about signs from index re-orderings. One 
has to note that certain problems are specific to situations with quite many 
indices, so that simple test cases frequently just do not exist. The problem 
may well be that even the appropriate questions have not yet been found that 
have to be asked, and answered, to make the validation phase more bearable. 
Addressing this problem presumably should be the most reasonable next step 
before changing the definition of normal forms such that Garnir symmetries are 
properly taken care of. Then, the issue of a systematic determination of [E, F] 
commutators has to be addressed before the Eiq sigma-model can be extended 
to higher levels. As it turned out that the task of actually implementing the 
relevant algorithms is highly prone to sign and similar errors, and as debugging 
is complicated by the fact that it is virtually impossible to develop an intuition 
for such terms, a reference implementation in LISP that also can be used to 
redo the calculations presented here within reasonable calculation time will be 
made publicly available in the next release of |1(J|. 
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A Algorithms 



In this appendix we list algorithms that are used to address various parts of the 
problem. 

A.l Efficient index (anti-)symmetrization for tensors of 
high rank 

The naive translation of the symmetrization/anti-symmetrization prescriptions 
for Young tableau projections to a computer algorithm reaches its limits quite 
early, often already with problems that - using some thought - easily can be 
done by hand. The primary underlying reason is that a naive execution of a 
task like e.g. the anti-symmetrization over nine indices generates 9! = 362 880 
different individual copies of a term with indices re-distributed in all possible 
ways which then usually have to be classified and subjected to further reduc- 
tions. Nevertheless, quite many research problems where such manipulations 
are necessary are of such a structure that, given an unlimited amount of time 
and energy, the analysis could be taken to almost arbitrary depth. Not sel- 
dom, there is a considerable gap between the point where one can claim to have 
obtained a sufficiently deep understanding of the mathematical structure that 
going any further is not expected to produce valuable additional insights and 
the point where it becomes unfeasible to do calculations by hand. 

This raises the question how to formalize those ideas that make the simple 
cases feasible by hand in such a way that they can be executed as effectively 
by a machine, only much faster and much more reliable than by a human. The 
key idea here is to make as much use of overlapping symmetrizers and anti- 
symmetrizers as possible. 

A. 1.1 Data Structures 

While one may want to consider terms formed out of general products of tensors 
with mixed symmetry, we limit our discussion to the framework of Lie algebras, 
where the only product of terms is the commutator. Conceptually, this is also 
the most interesting the issue of normal forms of terms will be more 

involved here than with simple tensor products. Additionally, we restrict our- 
selves to tensors with contravariant indices belonging to one SL{n) algebra only. 
This last restriction will eventually have to be softened, at the very least to al- 
low indices to be either all contravariant or all covariant, in order to represent 
commutators of the [E, F] type. The case of more general base algebras than 
SL{n) may also be studied, but the price that has to be paid is that the relation 
between irreducible representations and corresponding tableaux gets much more 
involved, hence this is not considered here. 

Tensors The most natural way to represent a tensor T is as a tuple T = 
(r, u), 7, 4>), where r is the type of the tensor, u> is some abstract weight that will 
be used to normalize commutators, 7 the shape of the Young tableau, and <f) a 
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vector of indices. On the pragmatic level, r will also include information about 
a printable tensor name. We furthermore assume a strict order u> < ui' as well as 
addition to be defined on tensor weights, which will be used for sorting tensors, 
e.g. to determine normal forms for sums. The tableau shape 7 must provide 
information about where to place an index from the linearized presentation <j) in 
the tableau. This can be achieved e.g. by just recording the lengths of the rows 
of the tableau, and using the convention to record indices in normal reading 
order. Furthermore, 7 must provide information that allows one to quickly 
extract all the indices belonging to a given column from </>. When building 
terms from tensors, it is useful not to put the tensor T into the term as it is, but 
rather a tuple T» = (T, n) , where the tensor is supplemented by an extra bit of 
information describing whether reduction to normal form has been applied to 
this quantity or not. Note that data representations in a real implementations 
may not strictly adhere to these lines, as there might be practical reasons to 
e.g. include normalization information directly in the tensor data structure. 

Terms From tensors, we can recursively form sums and commutators. This 
is formalized by the definition of a term S as a tuple S — (s, £), where s is a 
(sorted) vector of summands {K, a) with K either a tensor or a commutator, 
and a a coefficient. We impose the restriction that all summands must be of the 
same weight. The weight of a commutator is just the sum of the weights of its 
summands, and the weight of a term is just the weight of any of its summands. 
£ describes all the non-overlapping symmetrizations or anti-symmetrizations 
applied to the indices occurring inside this term. As index-symmetrization is an 
operation that is performed by multiplying a 5^ " tensor onto an expression, 
this is most naturally represented on the level of terms, not individual tensors. 
In particular, symmetrizations will frequently stretch over indices belonging to 
two different sides of a commutator. £ is a list of pairs a = (p, v) of a type 
bit discerning symmetrization from anti-symmetrization and a lexicographically 
ordered vector of index names. £ itself is ordered lexicographically with respect 
to the index name vectors. 

As we may occasionally have to deal with the problem to extract a sign 
factor from a term, e.g. when bringing a commutator to normal form, which 
would have to be propagated through the 'outer context' of a subterm that is 
embedded in a term for normalization purposes, it can be convenient to be able 
to define normal form modulo a constant factor, hence we consider augmented 
terms = (5, /, n) which arc tuples of a term S with coefficient of the leading 
summand being 1, an overall factor /, and a bit denoting whether S and all its 
summands are guaranteed to be in normal form. 

Commutators A commutator C is just a pair (£, r) of a left and a right 
side, which both are augmented terms £ — Si*,r — 5*2*. Again, we augment 
this structure by an additional bit, C* = (C, n), that provides information 
whether both the left and right entry are guaranteed to be in normal form, and 
furthermore whether the left side is smaller than the right side with respect to 
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the order on terms, to be described next. 
A. 1.2 Normal Forms 

Consequent use of augmented tensors, commutators, and terms allows to do 
many symbolic manipulations without unnecessary normalization steps for sub- 
terms that already are in normal form. Just as the definitions of commutators 
and terms have to recursively refer to one another, so do the definitions of their 
respective normal forms. We define a term to be in normal form, if its sum- 
mands are in normal form and sorted in ascending order with respect to the 
order that treats a commutator as smaller than any tensor, compares commu- 
tators lexicographically, and orders tensors first by their type r, then by their 
index distribution (which is assumed to be normalized), and read for the pur- 
pose of ordering in Kanji reading order (column by column right to left, every 
column up to down). A commutator is in normal form if the summands of the 
left side are smaller in lexicographical comparison to the summands on the right 
side. For index symmetrizations, normal form is defined by first ordering the 
indices in every symmetrization block lexically, then ordering symmetrizations 
lexicographically by their indices (not symmetrizing/anti-symmetrizing type). 

Commutator and tensor normalization is complicated by the fact that these 
may be subject to index symmetrizations that do not match their own structure 
(i.e. indices to be antisymmetrized over may be spread out over both sides of 
a commutator, or multiple columns of a tableau). This is resolved as follows: 
when recursively traversing the tree of commutators and terms, and mapping it 
bottom-up to a tree in normal form, an environmental parameter that denotes 
the presently active index symmetrizations is set whenever recursion proceeds 
from a term with extra symmetries to its summands. In some programming 
languages like Perl or LISP, this is facilitated by the use of local variables, 
which, when set in the caller are visible to the callee and shadowed by inner 
local definitions. This can prove to be useful in order to unclutter code from 
extra function arguments that perhaps should be kept implicit. 

Tensors are normalized by first normalizing all extra index symmetrizations 
that act on it to those indices that occur in the tensor, then applying them se- 
quentially to bring as many lexically small indices to the rightmost column, then 
to the second rightmost, etc. Then, tensor columns are sorted lexicographically, 
and extra indices protruding in the first column over the second, which are only 
subject to symmetrization from the Young tableau projector, are brought to 
inverse lexical order (which will again move the lexically smallest index to the 
rightmost column). Clearly, great care has to be taken with signs introduced 
by these symmetrizations. 

With the normalization of commutators, there are two aspects: first, which 
of the terms to place left, and second, how to distribute indices between them 
should an index symmetrization from the environment cover indices on both 
sides? Unfortunately, these are somewhat interlocked, and while this issue may 
be resolved properly, this would mean both unnecessary code complexity and 
extra calculation time for many applications. Hence, we use a somewhat sim- 
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plistic approach to this and first compare the weights of the terms on both sides. 
Should these turn out to be equal, we do a lexicographical comparison to the 
summand vectors, with the same order that is used to order summands. This is 
then used to decide which term to put to the left (the one with larger weight, or 
the second in lexicographical comparison for equal weights) in the commutator. 
Then, we make use of the presently relevant symmetrizations one after another 
to rename the indices in the commutator in such a way that each one tries to 
bring as many lexically small indices to the left as possible. 

One should keep in mind that with this definition of normal forms, applying 
a Young symmctrization to a term containing a single tensor that is in normal 
form will in the general case produce a term that contains more than one tensor 
in normal form, as in: 
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A. 1.3 Symmetrizations 
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a(abcdefghi) ■ a(jkl) 
(A.34) 



The basic operation is to apply an extra index symmetrization to a term which 
typically will overlap with some of those already present. To give a specific 
example: 
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(A.35) 

Taking a term S that carries N symmetrizations <jj = (pj,vj),j e {1, . . . N} 7 
where we take the Vj to cover all of the indices of S (i.e. we include anti- 
symmetrizations over single indices for uniform treatment), it is comparatively 
easy to introduce a new symmetrization a' that completely covers a subset of the 
existing ones: the new index symmetries of the re-symmetrized tensor are just 
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the pre-existing ones plus a' minus those absorbed into a', brought to normal 
form. One only has to recursively traverse the term to bring it back to normal 
form with the new index relabeling symmetries in place. 

For the case that a new symmetrization a' will only overlap some of the 
pre-existing symmetrizations without covering them completely, a reasonable 
approach is to first split those overlapping symmetrizations one by one, and 
then proceed as above. 

Splitting symmetrizers is such an important and ubiquitous function that 
a symbolic algebra library should export it to the user. While one normally 
would want all such functions to return normalized terms, it turns out that 
there are a few situations where such normalization is not immediately helpful, 
but computationally expensive, hence it should be possible to disable the final 
normalization step with an appropriate parameter. 

Besides the obvious term to be split, the further arguments to the index 
splitting function are the previously existing symmetry a — (p, v) and a (sorted) 
vector of sub-indices v' . Re-distribution of indices in all possible ways modulo 
permutation of those n\ indices that are in v' and permutation of those n 2 
indices that are in v and not in v' generates - with our normalization conventions 



extra coefficients ±l/c, and - in the latter case - if the new index re-distribution 
is an even permutation of the former one or not. 

Assuming the availability of an iteration construct (to be discussed later) 
that for a given pair of natural numbers (fci,fc 2 ) implements systematic pro- 
cessing of all different pairs of vectors (101,102) of respective lengths fci, fc 2 for 
which the concatenation W1W2 is a permutation of the vector of natural num- 
bers [0, . . . , ki + fc 2 — 1] and for which W\,W2 are each sorted in ascending order, 
index permutation splitting works in detail as follows: 



h be a modifiable (hash) table mapping summands to coefficients, 
£' be the list of symmetrizations, 

with a removed and (p, v'), (p, v — v') added instead 
p a be a vector of m natural numbers 

such that p® is the position where the j-th index of v' 

appears in v, 
p b be a vector of n 2 natural numbers 

such that p b is the position of the j-th index of v 

that does not appear in v' , 
a be the concatenation p a p b , 

and S' be a copy of the term S with its associated symmetrizations 




contributions obtained by suitable index renamings, with 



let 



replaced^ by £' 



then 



For all c ordered choices (wi,w 2 ) of vectors of natural numbers 
of lengths #W\ = m, #w 2 = n 2 , with w = jiji!o 2 , 
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let 

v* be the vector with entries v* — v a . , 

7 be the permutation factor, ±l/c, where the negative sign 
only occurs for antisymmetrizations and 
when v* is an odd permutation of v, 

S v * be the term S' where indices are relabeled^**) 
such that Vj \— > Vj, 

then 

add all the summands in S v * multiplied with 7 

(also taking the overall factor / from S v * into account) 
to the result table h. 

then 

Generate a result term from h, 

and normalize unless expressly wanted otherwise. 

Remarks: 

(*) This must happen here since we have to avoid working on terms with 
inappropriate symmetrizations attached to them in subsequent steps! 

(**) One can obviously omit performing the relabeling of symmetrizers here. 

A. 2 Generating all bare level £ commutator structures 

Given a set of n Lie algebra generators X\, X 2 , . . . , X n , we want to generate all 
recursive commutator structures from them that are unique up to application 
of [A, B] = — [B, A]. For example, there is only one such structure at level 1 = 2 
(LY X ,X 2 ]), while there arc three at I = 3 ([[X lt X 2 ], X 3 ] ,[[X U X 3 ] , X 2 ] ,[[X 2 , X 3 ] , X 1 
and 15 at £ = 4 (three of type [[X a , X b ] , [X c , X d }] and twelve of type [[[X, X],X], X]) 
We inductively define a strict order on commutator structures P, Q: 



P<Q 



P is one of the X n , and Q is X m , and n < to, or 
P is a commutator, and Q is one of the X n , or 
P and Q are commutators [P a ,Pt,] and [Q a ,Qb] 

and the vector (P a ,Pb) is lexicographically 

smaller than (Q a , Qb) 



(A.36) 



We consider a commutator structure to be in normal form if all the commu- 
tator structures contained in it arc in normal form, and it is either one of the 
X n or a commutator [P, Q] where P < Q. 

Given that all commutator structures in normal form of level up to £ = N 
are known, we can form all commutator structures of level up to I = N + 1 by 



taking all commutators [P^—j, Qi=jv+i-j] °f au lcvel-j commutator structures 
P' with all commutator structures Q' of complementary level N + 1 — j smaller 
than j, where P',Q' arc formed from P, Q by suitable relabeling of indices. 
The P, Q — ► P',Q' index relabeling is done in two steps: first, one offsets all 
indices on the X a in Q by j, then one performs all possible ^ N j~ 1 ) index 
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permutations that preserve the relative order among the indices in the range 
1 . . . j and j + 1 . . . N + 1. For N + 1 even, there are extra commutators with 
P and Q of the same level, but one has to take care to take only those after 
relabeling where [P 1 , Q'] is in normal form. 

The number of different commutator structures grows as follows with the 
level: 



Level 


^Structures 


1 


1 


2 


1 


3 


3 


4 


15 


5 


105 


6 


945 


7 


10 395 


8 


135135 


9 


2 027 025 



(A.37) 



A. 3 Generating all bare level £ Jacobi-type relations 

The commutator structures generated by the previous algorithm are not in- 
dependent. While they honor the fact that [A, B] and [-8,^4] are dependent, 
they do not yet take into account the Jacobi identity [[A, B],C] + [[£>, C], A] + 
[[C, A],B] = 0. We want to systematically generate all such identities between 
three commutator structures that can be obtained by applying such a 'Jacobi 
rotation' to any double commutator. The list of identities generated in this step 
may not be algebraically independent, but will be complete and not contain 
duplicates. 

For an efficient computer algorithm, it is important to make use of the 
distinction between 'the same' (value in the computer's memory) and 'equal' 
(with respect to comparison). While it is always easy to find out if two values 
are the same, finding out of they are equal may be costly for recursively built 
trees, hence the algorithm will contain canonicalization steps, where a freshly 
generated commutator structure is mapped to the one reference value in memory 
to which it is equal. We assume our memory representations to be such that 
two Xi are equal iff they are the same value, while two commutators that were 
generated (hence, memory-allocated) independently at different steps in the 
algorithm will be assumed non-equal. 

The input of the algorithm to generate all bare level t Jacobi relations is a 
list of all level I commutator structures, as generated by the previous algorithm. 
We in particular assume that no two commutator structures in this list share 
any memory representations for commutators. (Should the implementation of 
the previous algorithm not give that guarantee, we first perform a recursive copy 
on all structures to make every commutator unique.) 

In a first step, we create a dictionary C that maps (with respect to being the 
same, not equality) commutator structures to pairs of an unique tag (a number) 
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and the corresponding level. This is filled by recursively walking through all the 
input commutator structures. This means that for level 3, this table would hold 
3 • 2 entries: the three top-level commutators, and the three inner commutators. 

Furthermore, we create a dictionary M that maps (with respect to being 
equal) level i commutator structures to their canonical memory representation, 
which is just the input data. (For level 3, there would be three entries in 
this table.) A third dictionary T is used to remember (with respect to being 
the equal) pairs of tags of outer/inner commutators appearing somewhere the 
canonical memory representation of commutator structures that already have 
been used in some Jacobi relation. 

We then process the input commutator structures one after another, recur- 
sively walking through the tree of commutators, where at every place where 
we encounter a double commutator of either the form or [x, [y, z\], we 

obtain the tags of these commutators from C and check if these two generators 
were already used previously in the generation of a Jacobi relation by check- 
ing whether there already is a corresponding entry in the tag-pair dictionary 
T. If so, we do nothing at this position and just recurse further down. If not, 
we in addition make a corresponding new entry in T and form the other two 
sub-structures obtained by cyclically rotating x, y, z. These are both brought to 
normal form, taking care of signs. Note that the subtree depth that is stored in 
conjunction with the tag in C is useful to speed up determination of this normal 
form. 

The problem now is to first form the complete modified level i commu- 
tator structures for these, and second bring everything to canonical memory 
representation. A particularly useful technique to achieve this that is directly 
available in functional programming languages (and has to be emulated in other 
languages) is continuation coding: when recursively traversing the commutator 
tree A down from the root, we hand over two functions as arguments in the 
recursive call to the function processing subtrees S: the first is a buildup func- 
tion fis that, when called with a a commutator structure tree S', will generate 
a commutator structure that equals A with the subtree S replaced by S". The 
second is a selector function as that, when called with a commutator structure 
tree X will return (if possible) the subtree S of X that is situated in the same 
position relative to the root of X as is £ relative to A. In the recursive call 
to the left commutator structure subtree L of a commutator structure 5, the 
function arguments /3l,ctl are: 

0l = (£' ~ /?s(pair(L', right(S)))) 

o~ l = (X^Mt(a s (X))) [A - 6S> 

and likewise for the call to the right commutator structure subtree. 

First, the buildup function (3$ is used twice to generate the two other Jacobi- 
rotated level i normalized commutator structures. These are then mapped to 
their canonical memory representation by use of the dictionary M. The selector 
function as now allows to locate the canonical memory representations of the 
rotated outer and inner commutator. These are looked up in C to form their 
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tag pairs, which are then entered in T to ensure that the relation that just has 
been found is not generated again (twice, even). 

The number of relations grows as follows with the level: 



Level 


#Jacobi relations 


1 





2 





3 


1 


4 


13 


5 


150 


6 


1818 


7 


25 221 


8 


400 260 



(A.39) 



For level I = 4, we get explicitly: 



= - 




,x 3 ] 


,X 4 ]- 


'[Xi,X 2 ] 


[x 3 ,x 4 \]- 


[[[Xi 


x 2 ] 


,x 4 ] 


,x 3 ] 


= - 


h[[[*i>*4 


,x 2 


,x 3 ] - 


[Xi , X 4 \ 


[x 2 ,x 3 \]- 


[[[Xi 


x 4 ] 


,x 3 ] 


,x 2 ] 


= - 




,x 2 


,X 4 ] - 


[Xi,X 3 ] 


[x 2 ,x 4 ]\- 


[[[Xi 


x 3 ] 


,x 4 ] 


,x 2 ] 


= - 


^[[[x 3 ,x 4 


,x 2 


,Xi] - 


[[X 2 ,X 4 


7 X 3 ], Xl] + 


[[[X2 


x 3 ] 


,x 4 ] 


,Xi] 


= - 


^[[[X 3 ,X 4 


,Xi 


,x 2 ] - 


'[[Xi,X 4 


7 X 3 ] 7 X 2 ] + 


[[[Xi 


x 3 ] 


,x 4 ] 


,x 2 ] 


= - 


^[[[x 2 ,x 4 


,Xi 


,x 3 ] - 


'[[Xi,X 4 


7 X 2 ], X 3 ] + 


[[[Xi 


x 2 ] 


,x 4 ] 


,x 3 ] 


= - 




,Xi] 


,X 4 ] - 


'[[Xi,X 3 


, X 2 ] , X 4 ] + 


[[[Xi 


x 2 ] 


,x 3 ] 


,x 4 ] 


= - 


H[Xi,X 4 ] 


[x 2 


x 3 }] + 


\[[Xi,X 4 


7 X 3 ] , X 2 ] — 


[[[Xi 


x 4 ] 


,x 2 ] 


,x 3 ] 


= - 


H[Xi,x 4 ] 


[x 2 


x 3 }}- 


[[X 2 ,X 3 


7 X 4 ],Xi] + 


[[[X2 


x 3 ] 


,Xi] 


,x 4 ] 


= - 


H[XuX 3 ] 


[x 2 


X 4 }} + 


'[[Xi,X 3 


7 X 4 ],X 2 ] — 


[[[Xi 


x 3 ] 


,x 2 ] 


,x 4 ] 


= - 


H[Xi,x 3 } 


[x 2 


Xi}} - 


[[X 2 ,X 4 


7 X 3 ),Xi] + 


[[[X2 


x 4 ] 


,Xi] 


,x 3 ] 


= - 


H{XuX 2 ] 


[x 3 


X 4 }} + 


'[[Xi,X 2 


,X 4 },X 3 ] — 


[[[Xi 


x 2 ] 


,x 3 ] 


,x 4 ] 


= - 


V[[Xx,X 2 ] 


[x 3 


Xi}} - 


[[x 3 , x 4 


7 X 2 ],X{\ + 


[[[x 3 


x 4 ] 


,Xi] 


,x 2 ] 



These equations then have to be instantiated with X\ = Ef bc , X 2 = Ef e ^, 
etc., before definitions of known commutators are applied to obtain relations 
between commutators and tensors at the given level. Then, all different ways 
how to distribute indices from the groups abc, def, etc. to the columns of a 
given Young tableau have to be found (modulo exchange of groups). These 
then have to be brought to the corresponding tensor's anti-symmetries in index 
columns. 
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